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1. INTRODUCTION 

During recent years, a lot of effort has been devoted to studying the space-time 
evolution of ultrarelativistic heavy ion collisions. Especially, from the point of view 
of quark-gluon plasma (QGP) formation, it is important to understand energy deposi- 
tion and particle production mechanims in the central rapidity region of the collisions. 
Qualitatively these mechanisms can be divided into two categories: semihard and soft. 

Semihard particle production is based on perturbative QCD (pQCD). The semihard 
gluons, quarks and antiquarks are produced in parton-parton collisions and they have 
transverse momenta pt > Po ~ 1...3 GeV. These semihard QCD quanta are often called 
minijets [1] because even in pp{p) collisions they are not observable as individual jets 
below p T ~ 5 GeV [2]. However, even if minijets are not at all directly observable in 
nuclear collisions, they are expected to play a major role in the space-time evolution of 
central collisions at collider energies [3-7] Below the momentum transfer scale ~ 1 — 2 
GeV pQCD is not expected to be reliable anymore, and modeling for soft particle 
production from beam jet fragmentation is needed. In heavy ion collisions this can be 
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done e.g. in terms of strings as in hadron-hadron collisions [8] , or in terms of a strong 
background color field decaying into particles through the Schwinger mechanism [9]. 

In the space-time evolution of an ultrarelativistic heavy ion collision the semihard 
and soft particle production mechanisms take place at different time scales. Semihard 
QCD quanta in the central rapidity region are typically formed at very early times, at 
r ~ 1/pt ^ 1/po ~ 0.1 fm, whereas the characteristic time scale for the soft processes 
(like the decay time of the background color field), is r ~ 1 fm. In this sense, the minijet 
system is formed first and it will then serve as initial condition for the further evolution 
of the system. An example of this kind of modeling can be found e.g. in ref. [10]. 

In the central rapidity region the relative strength of the semihard and soft particle 
production mechanisms depends strongly on the cms energy of the heavy ion collision. 
The semihard part is expected to play a more important role with increasing energy, 
so that it is negligible at energies y/s ^, 20 AGeV, becomes relevant at the Brookhaven 
Relativistic Heavy Ion Collider (RHIC) energies y^s ~ 200 AGeV, and finally dominates 
at the CERN Large Hadron Collider (LHC) energies ^/s ~ 5.5ATeV [3,4,6,11]. This 
conclusion is based on the idea of multiple minijet production as the cause of the rapidly 
rising inelastic and total cross sections of pp and pp collisions at energies > 200 GeV. 
A detailed discussion and further references can be found e.g. in ref [11]. 

Since the first ideas of the importance of the semihard processes in ultrarelativistic 
heavy ion collisions [3,5], a lot of inspiring work has been done on creating event gener- 
ators like HIJING [6] and Parton Cascade Model (PCM) [7]. Naturally, the techniques 
of treating semihard processes in event generators like PYTHIA [12] for hadron-hadron 
collisions has been of great help in building up heavy ion event generators. In spite of 
the differences in initialization and in final state interactions, the semihard processes 
lie in the heart of generating events with e.g. HIJING and PCM at heavy ion collider 
energies. Therefore, to have a better control over the results from simulations, it is im- 
portant to understand the basic processes as well as possible, and especially to study the 
uncertainties, like nuclear and hadronic shadowing of the parton distributions, higher 
order contributions to the cross sections, and the cut off scale po used to give the division 
between semihard and soft physics. 

In this talk, I will concentrate on the semihard parton production only. I will first 
estimate the initial conditions for early QGP formation. The main emphasis will be in 
consequences of the small-x enhancement of the parton distributions observed in the 
HERA measurements, as studied in ref [13]. After considering the uncertainties in the 
pQCD calculation, I will discuss the evolution of the formed QGP with minijet initial 
conditions, and finally I will give an estimate for the lower limit of multiplicities for 
LHC nuclear collisions. 
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2. FORMATION OF QGP 



2.1. Impact of the HERA measurements on minijet and Et production in 
AA collisions 

The new deep inelastic ep-data from HERA [14] exhibit an increase in the structure 
function F^ix^Q ) at small x relative to the behavior implied by earlier data. Distri- 
bution function analyses using the new data have been carried out. These analyses, in 
particular, constrain the gluon distribution function, which is not directly measured. 
The estimates to be presented in this talk are computed with parton distributions DO', 
D-' and H by Martin, Roberts and Stirling (MRS) [15,16]. The corresponding gluon dis- 
tributions for these sets at Q = 2 GeV are shown in Fig. 1. At small x the distributions 
shown behave as x~ s with 5 = 0, —0.3 and —0.5, for DO', H and D-', correspondingly. 
In the global analysis by MRS, the set H follows from the best fit to the HERA data 
[16]. 




Fig. 1. Gluon distribution functions MRS DO', D-, and H [15,16] at the scale Q = 2 
GeV. The old gluon distribution of Duke-Owens set 1 (DOl) [17] is shown for compar- 
ison. For the vertical lines, see the text. 



For minijet production at high energies the consequences of the observed small- £ 
rise are obvious. Minijets falling into the mid-rapidity y = dominantly have transverse 
momenta px ~ Po ~ 2 GeV, and they mainly come from momentum fractions x ~ 
2p /yfs. In Fig. 1, the vertical lines illustrate the typical value of x for LHC and 
RHIC with y/s = 200 GeV and y/s = 5500 GeV, correspondingly. Difference between 
the MRS distributions shown is relatively small at RHIC, whereas at LHC there is a 
sizable difference between the sets, and this is reflected into the minijet cross sections 
discussed below. 
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When considering the number of produced semihard gluons, quarks and antiquarks, 
both in hadronic and nuclear collisions, the key quantity is the integrated (mini)jet cross 
section 

where the inclusive differential jet cross section is given by 

J ° 1 = K ^ x 1 f i/p (x 1 ,p^)x 2 f j / p (x2,PT)-^ ■ (2) 



dp\d yi dy 2 ^ /P "' v ~"" J 'dt 



ijk 



The momentum transfer in the process r/' — > kl is px 7 the fractional momenta of the 
partons i and j are #1,2, and yi^ are the rapities of the partons k and /. The calculation 
is only done to leading order (0(or s )). The higher order terms are in principle available 
from the calculation of S.D. Ellis-Kunszt-Soper (EKS) [18] but here these are simulated 
simply by a constant factor K ~ 2, which seems to be an acceptable approximation 
with a scale choice Q = pr [19]. 

In the eikonal formulation (see e.g. [11] and references therein) we can interpret 
2o"jet/c"hf as the average number of semihard partons with px > Po produced in an 
inelastic pp collision. This ratio may become much larger than 2, (see e.g. a- iet shown 
in fig. 1 in [13]) especially at LHC energies with small- a; enhanced gluon distributions, 
which indicates that one is effectively including minijet pairs not only from multiple 
(independent) semihard collisions but also contributions from the initial partonic flux, 
i.e. from the gluon ladders leading to the virtuality Q ~ pr [20]. However, in the central 
rapidity, the pairs should dominantly come from the semihard collisions. 

In the first approximation, where all nuclear modifications to the parton densities 
fi are neglected, the average number of minijets produced in an AA collision with an 
impact parameter b is obtained from 

N AA (h, vG,p ) = T AA (b) 2(7 je t(v^,Po) (3) 

where T AA (b) is the standard nuclear overlap function. Here we will need only T AA (0) « 
A 2 j \tx R a ) ~ 32/mb for A = Pb. Notice that there are two reasons for large N AA for 
heavy nuclei: the overlap function scales as ~ A 4 / 3 and the perturbative cross section 
Ojet becomes large at high energies. 

The results for Oj e t integrated over all rapidities can be found in [13], but what is of 
greater interest for us here, are the numbers of minijets falling into the central rapidity 
unit \y\ < 0.5. By making this cut in Eq. (1) (as will be done in Eq. (4)) and by careful 
bookkeeping of the flavors which fall into this rapidity window, we find out an estimate 
for the chemical composition of the produced minijets at LHC and at RHIC. This is 
shown in Table 1. 

We make two observations: First, as expected, the small-x rise of the parton dis- 
tributions does not affect much the numbers for RHIC, but inreases considerably the 
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Table. 1. Values of Naa(0, y/s,po) calculated for p = 2 GeV, y/s = 200 GeV, 
5500 GeV and for the three parton distributions considered in the text. The upper 
numbers are for LHC and the lower ones for RHIC. The results with Duke-Owens set 
1 distributions (DOl) are shown for comparison. An overall factor K = 2 but no shad- 
owing is included. Both the total as well as the contribution from gluons, quarks and 
antiquarks is given. 



numbers for LHC. Secondly, the formed minijet plasma is obviously strongly gluon dom- 
inated: about 80 % of the produced partons are gluons. At this point, one should bear 
in mind that the figures above are only the lowest order estimates, and the number of 
partons is well defined. However, when considering higher order terms in the jet cross 
section, one encounters collinear singularities in the case of three-jet final states, and, 
without a definition of a jet size R in the (y, (/>)-plane it is not clear when the partons 
should be counted as separate jets and when as one jet [18]. In general, radiative effects 
would tend to increase the gluon dominance of the formed minijet plasma. 

Let us then consider perturbatively the transverse energy in the central rapidity 
region of AA collisions. This can be estimated (see [4] for details) in terms of the first 
Pt moment 

/ /— \ 1 7-i \ f 1 d(T i e t / \ / . \ 

0jet(Vs,Po)(-ET) = / dp T dy 1 dy 2 — - 1 — (ei + e 2 )pr, (4) 

Jpo 27r 2 dp T dy 1 dy 2 

where the rapidity cut is made by the functions e^: ejt = £(PT,Vk) = 1 if the parton k 
falls inside the region of our interest and =0 otherwise. To illustrate the p dependence, 
crjet(-E'T) is plotted as a function of po in Fig. 2. Note again that no shadowing is 
included in the calculation. The composition into contributions from gluons, quarks 
and antiquarks with p = 2 GeV is shown in Table 2. The estimate for the total initial 
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Et in one unit of y near y = is then 



^ &Pb (0)(v^,Po, \y\ < 0.5) = T Pb p b (0)a iet ( v G,po)(E T } 

= 3200 (387) GeV DOl 

= 5440 (563) GeV DO' 

= 10300 (547) GeV H 

= 17600 (538) GeV D-' 

where the numbers are for LHC (RHIC), correspondingly. 
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Fig. 2. The first pt moment crj e t(v / i, Po)(Et) from Eq. (4) with the sets of distribution 
functions marked on the figure and for y/s = 200 and 5500 GeV [13]. 



Again, the increase in E T associated with the small-x increase is significant for 
LHC. However, this Et is the final Et only if no further interactions take place in the 
system. The largest decrease in Et is obtained if the system thermalizes and an ideal 
hydro flow is generated. Another interesting point is that in spite of the apparently 
large values of Et in Eq. (5), at most only about 1.5% of the total energy is stopped 
perturbatively. 

Since the longitudinal size of the comoving volume is Az = rAy and Ti ~ 1/po = 
0.1 fm the obtainable energy densities can be estimated with the Bjorken formula [21] 
as e{ji) = (E^(\y\ < 0.5))po/nR 2 A . Already without the smalkc enhancement of the 
parton distributions, i.e. for distributions like MRSD0', the energy densities e(r;) be- 
come large as compared to the critical density e c ~ 1...3 GeVfm -3 associated with 
the deconfinement phase transition. For MRSD0' we get e(rj) ~ 410 (43) GeVfm -3 
for LHC (RHIC). An interesting observation is that for distributions like MRSD0' at 
high enough energies the initial energy density does not depend strongly of po- This is 
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Table. 2. As Table 1, but for a- ]et (j>o)(E T ) in units of mbGeV with \y\ < 0.5, y/s = 
200 GeV, 5500 GeV and p = 2 GeV. 

because of the scalings: Oj e t(po) ~ Po 2 ■> a )et(Po)(ET) ~ Po \ an d the po-dependence in 
eiji) is canceled out. With the small- a; enhancement of the parton distributions, the 
results depend more on p 0: as illustrated in Fig. 2. 

About the approximations done, one should note that the transit times of the 
colliding nuclei were neglected. In AA collisions with A ~ 200 at RHIC the transit time 
is about 0.2 fm, which will reduce the estimate for the maximum energy density by a 
factor ~ 2 [10]. For LHC, the transit time is negligible as compared to 1/po = 0.1 fm, 
but nuclear shadowing may reduce the estimate at least by a similar factor (see e.g. 
[34]). 

2.2. Chemical equilibrium at r = 0.1 fm? 

The lowest order pQCD results of [13] described in the previous chapter show that 
the minijet-gluons dominate the early QGP. Let us now study how far from chemical 
equilibrium we actually are initially at r ~ 0.1 fm. 

For an ideal gas of massless gluons, quarks and antiquarks, we have the following 
ratios at zero chemical potential: 

32 



n g 



= ^< 0.9, 

ideal 9N f € 1 + ^ 



ideal 21AT, 



< 0.8, N f > 2 (6) 



n q + Ug 

From the previous lowest order pQCD results, we can extract the ratios as 

e 



n g 



n q + riq 



-3.7 (RHIC) 5(LHC), 

pQCD e Q + C <7 



-3.5 (RHIC). ...5(LHC) (7) 

pQCD 

The conclusion is the same from both number and energy densities and for both RHIC 
and LHC: the pQCD ratios are far from the equilibrium ones. Many more quarks and 
antiquarks would be needed for complete chemical equilibrium. 
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Then, an interesting question is that if we consider gluons only, can they be in 
chemical equilibrium, i.e. thermalized at r ~ 0.1 fm? In this case it would mean that 
the gluon number density is close to the thermal density obtained by assuming that all 
the initial energy density, e g (Ti), is fully thermalized. Therefore, for LHC, we compare 
energy/gluon computed in two ways: 

i) take the previous pQCD results from tables 1 and 2 to get E Pb g Pb and Np bPb with 
po = 2GeV, \y\ < 0.5. 

ii) take the pQCD estimate for E^ b g Pb : convert this into equilibrium (p, g = 0) energy 
density e g (ri) = 3a7i 4 , where a = 167r 2 /90, solve for T;, and get energy/gluon as 
e g( T i)/ n g( T i) ~ 2.77V For LHC, the results for different parton distributions are 

£/gluon = 4.1 GeV 2.7T, = 2.0 GeV DOl, 
= 3.6 GeV = 2.2 GeV DO', 

(8) 

= 3.2 GeV = 2.6 GeV H, 

= 2.9 GeV = 3.0 GeV D-'. 



For the distributions with xg(x) ~ const the thermalization of gluons requires a degra- 
dation of their average energy by collisions. For RHIC this is also the case [22]. With 
the HERA parton distributions the secondary interactions are only needed for changing 
the directions of momenta to make the distribution uniform. The conclusion is that the 
assumption of fast initial thermalization seems reasonable at LHC when new HERA 
parton densities are used. 



2.3. Uncertainties: K, p and shadowing 

In the previous results for the jet cross sections there are three main uncertainties: 
the K-factor, the momentum cut-off po and shadowing of the parton distributions. 

The calculations were done only to the lowest order (LO) pQCD, and a factor 
K = 2 was used together with a scale choice Q = pr to simulate the contributions from 
higher orders of a s . Here one should note that a 'K- factor' always depends on the scale 
choice; it can even be removed by choosing the scale differently. The Born level jet cross 
section does not depend on the jet size R in the (y, </>)-plane, while the next-to-leading 
order (NLO) calculation [18] does, and the experimentally measured jet cross sections 
do. Therefore, the X-factor always depends on the jet size, regardless of whether K is 
defined by comparing the LO results to the inclusive jet data [23] or to the LO+NLO 
calculation. Naturally, K also depends on the parton distributions and may also depend 
on the cms energy. 

Although I am not discussing the inclusive jet production here but, rather, consid- 
ering more semi-inclusive quantities, (i.e. all minijets and their transverse energy in the 
central rapidity unit), I show in Fig. 3 the inclusive 1-jet distribution from [19] in pp 
collisions at y/s = 200 GeV with Q = pr and R = 1. This shows that with this scale 
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Fig. 3. The inclusive one-jet cross section da/dprdy vs. pr at y = in collisions 
at y/s = 200 GeV, as predicted by the NLO (solid) and LO (dashed) calculations with 
Q = px- The MRSD-' parton distributions were used. For the NLO calculation R = 1.0. 
The figure is from ref. [19], and the calculation was done by using the program by EKS 
[18]. 



choice, it is a good approximation to have constant K ~ 2 to simulate the NLO terms 
in the jet cross sections. 

When the NLO contributions to the jet cross sections are included, the dependence 
on the scale choice should be weaker. This has been nicely demonstrated in [18] (see 
also [19]), and it works well for jets with pt > 10 GeV. However, when approaching the 
minijet region px ~ Po ~ 2 GeV, the convergence of the perturbation series becomes 
worse, and one would eventually need even higher order terms in the jet cross sections 
in order to improve the determination of the scale choice. This does not mean that the 
LO minijet results would be incorrect, it simply means that it may be very difficult to 
pin down the accuracy of the minijet calculation below a factor 2. As an example of 
this, we see in Fig. 3 that at least down to pr = 5 GeV the jet NLO calculation does 
not exhibit any striking behavior as compared to the LO result, even though at this 
low pt the NLO result depends on the scale choice practically as much as the LO does. 
However, it is clear that this problem needs further studies. 

In the approach considered here, the parameter p defines a division into semihard 
and soft physics. Therefore, po is a phenomenological parameter, and data on multi- 
plicities and cross sections are needed to determine p . In the eikonal approach for the 
pp(p) collisions, the cross section cr so ft for processes with pr < po is modeled in together 
with orjet(po)- Through the condition that experimentally measured pp(p) total cross 
sections are reproduced, the value of o"j et (po) then depends on the choice for (j so ft [1,11] . 
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In principle, in determination of po, one should also compare the final multiplicities as 
was done in [11]. Rather than doing a complete analysis like this, I have shown the 
results from [13] as a function of po- However, a reanalysis on determination of po with 
the new HERA parton distributions should be done, and the possible A dependence of 
Po should be studied in more detail. 

Related to the cut-off po, I note the following. Part of the higher order contributions 
come from 2 — > 2 + (n > 1) processes. At the tree level, these have been studied 
recently by E. Shuryak and Li Xiong [24], based on the matrix elements in [25]. In 
this treatment, one is faced with the jet resolution problem mentioned before, and with 
cut-off Sij > so ~ 4 GeV 2 for each pair ij, the higher order tree graphs seem to become 
dominant over the 2^2 processes at LHC energies [24]. Again, this shows that the 
actual number of gluons depends on the resolution, and is connected with color screening 
as well. Certainly, this interesting problem needs further studies, especially for LHC 
energies. 

Let us then come to the uncertainty due to modifications of the parton distributions. 
Since here we are mostly interested in the small-x phenomena, I will only consider 
shadowing. 

Shadowing in a proton means a depletion or a saturation in the gluon and sea 
quark distributions at x < 0.01. For the distributions behaving as xg ~ x~ s , at suffi- 
ciently small values of x and/or Q the partons start to overlap spatially, and a reduction 
in the parton density results due to the occurring recombinations. Corrections from the 
gluon fusions to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) scale evolu- 
tion [26] at small x have been calculated in the leading double logarithm approximation 
by Gribov-Levin-Ryskin and Mueller-Qiu (GLRMQ) [27]. In the semiclassical approxi- 
mation it has been shown that after correcting the evolution of the steeply rising gluon 
distributions ~ x~ 5 , they behave at asymptotically small x as xg ~ const [28]. At 
RHIC energies, the shadowing corrections in proton can obviously be neglected but for 
LHC energies they should be taken into account, at least with the new HERA parton 
distributions. 

There are two sources of gluon recombinations causing the shadowing in a proton: 
the two fusing gluon ladders, which couple four gluons to two gluons, can arise either 
from independent constituents of a proton or from the same one, as discussed in [27-29]. 
The HERA data for the structure function would seem to favor the first possibility; 
the latter source may be ruled out as too strong a fusion since no saturation in Fi is seen 
so far. However, since the sea quarks are emitted from the gluons they carry smaller x 
than the mother gluons, and the saturation of the quark distributions should show up 
at somewhat smaller x than that of gluons. It would be important to understand this 
better, since nuclear gluon shadowing arises perturbatively from the same sources as in 
the proton shadowing, with possibly different relative strengths [30] . 

Nuclear effects in the parton distributions in general have been widely studied dur- 
ing the recent years (see e.g. [31] for further references). However, a unified description 
through the whole x-range is very difficult to make, due to the different origin of the 
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effects at different ranges of x. From the point of view of phenomena discussed in this 
talk, shadowing as the small-x effect is the most important one. Notice also that an 
interesting theoretical idea for finding out the nuclear gluon densities at small x was 
discussed in this conference by R. Venugopalan [32]. 

Shadowing in a nucleus means a depletion of the nuclear parton densities rel- 
ative to the parton densities in a proton: f^(x,Q 2 ) < Af[(x,Q 2 ) at x<0.1. Unlike 
shadowing in proton, nuclear shadowing of the structure function F^ is experimentally 
clearly observed in deep inelastic scattering [33] , and it seems to be fairly independent 
of the scale Q. Although measured in fixed target experiments, nuclear shadowing gets 
a simple qualitative explanation in the infinite momentum frame in terms of overlapping 
partons. If the longitudinal wavelength of a parton, ~ 1/ (xp) , exceeds the contracted 
size of the nucleons (or the inter-nucleon distance) in the nucleus, ~ 2r n m n /p, shadow- 
ing should show up at x < 1 / (2m n r n ) ~ 0.1. According to the data this seems to be the 
case. Also, at sufficiently small x the partons from all nucleons at the same transverse 
location will overlap, and a saturation of f^/Aff is expected, at least when ff starts 
saturating (see e.g. [30]). 

Now, in the minijet calculation with the latest HERA parton distributions, in- 
clusion of shadowing is not a trivial matter. The nuclear shadowing effects are often 
simulated by simply multiplying the unshadowed parton distribution in a proton by a 
Q-independent ratio, extracted from the data [33] for F^/AFf . While this is presum- 
ably a rather good approximation for the sea quarks, it may not at all be the right 
thing to do with the gluons. Even when neglecting the shadowing in a proton - which 
should be a reasonable approximation at the RHIC energies - the ratio xgA/Axg p may 
be more strongly Q-dependent than that of the sea quarks. This was demonstrated in 
[34]. Besides the Q-dependence, one should also make sure that the valence quark and 
momentum sum rules are fulfilled for the nuclear parton distributions as well. 

At LHC energies the problem is more complicated since the gluon shadowing in a 
proton must not be neglected. Especially, with the rapidly rising gluon distributions, 
the gluon density in a nucleus at small x cannot be reliably approximated by simply 
multiplying by F^/AF^ since xg p should be shadowed as well. With the absence of 
a direct measurement for the gluon densities, the nuclear gluon shadowing can at the 
moment only be modeled in theoretically. By using the same perturbative mechanism to 
generate gluon fusions as described in the proton case above, nuclear gluon shadowing 
was studied in [30]. These results have not been implemented yet into the minijet 
calculation, but the procedure is clear: After extracting the nuclear parton distributions 
at some scale Q as in [30] and as in [34], one should then perform the full DGLAP- 
GLRMQ evolution in the case of nucleus. This work is in progress [35]. To get a first 
feeling of the order of magnitude of the shadowing effects, I refer to figs 5 and 7 of [34]. 
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3. EVOLUTION OF THE QGP 

Let us then briefly consider the early evolution of the QGP in a one-dimensional 
longitudinal scaling expansion with boost invariant initial conditions [21]. The initial 
conditions are given by the system of formed semihard gluons, quarks and antiquarks 
at t ~ 0.1 fm. As mentioned in the Introduction and as shown in [10], the soft processes 
are expected to be important for the evolution of the energy density e(r) up to RHIC 
energies but to be negligible at LHC energies (see e.g. fig. 2 of [10]). Since in this talk 
I am only considering the semihard contribution, let me now comment on the evolution 
of the QGP at LHC energies only. 

The QGP at r = 0.1 fm is strongly gluon dominated. Let us therefore simplify 
the discussion by neglecting the fermions. Furthermore, as demonstrated before, with 
the small- £ enhanced gluons it is reasonable to approximate the gluons to be fully 
thermalized at r = 0.1 fm. Then, the evolution of the energy density e g (r) of the gluon 
plasma is governed by the equation [21] 

UT T 

where P g is the gluonic pressure. If viscosity and Ohmic heating effects are taken into 
account, they would appear on the right-hand side of the equation above. These have 
been studied e.g. in [10] in connection with particle production from the decaying 
chromo-electric background field. However, for the discussion below these effects are 
not relevant. 

Now, there are two extreme cases for the evolution of the gluon plasma: 

i) free streaming: P g = =>- e g ~ r -1 

ii) ideal hydro: P g = e/3 => e g ~ r" 4 / 3 

The assumption of full thermalization of glue is relevant only for the latter case. If 
viscosity was included, it could at most reduce the PdV work to zero, so the first case is 
clearly an upper bound for the evolution of e g of a purely gluonic plasma (GP). In Fig. 
4 the evolution of e g is plotted as a function of r in the two extreme cases mentioned 
above. From the figure we see that if the deconfinement phase transition starts at 
e g ~ 1 — 2 GeV fm -3 , even the minimum life time of the GP at LHC is ~ 10 fm. 

3.1. Multiplicities at LHC - lower limit 

As a final topic, let me then estimate the multiplicities obtainable at LHC. The 
expected charged multiplicity is a crucial variable for the planning of experiments [36] 
and the various theoretical predictions lie roughly in the range 1500-8000. 

Certainly, the free streaming case above would give us an upper limit for the mul- 
tiplicity, ~ E PbPb (\y\ < 0.5)/0.5GeV. This would correspond to having no secondary 
interactions after the initial collision. In practice, this bound as such is of no realistic 
use, since the final state interactions are practically inevitable in this dense systems. 
The lower limit of multiplicity is obtained by considering only the gluons in the ideal 
hydro case (solid line in Fig. 4), where the system is fully thermalized at r = 0.1 fm 
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Fig. 4. The two extreme cases of the evolution of the energy density e g of the gluon 
plasma with initial conditions given by the pQCD calculation with the set MRSH [16]. 
Solid curve: ideal hydro flow with complete initial thermalization. Dashed curve: free 
streaming of the initial gluon system. 

and expands adiabatically thereafter. Then we only need to compute the initial entropy 
and convert that into the final state multiplicity. Pressure, energy density and entropy 
density are P g = aT 4 , e g = 3aT 4 , s g = 4aT 3 . Again, we estimate the longitudinal size 
of the comoving volume as Az = rAy and the energy density at the central rapidity 
region becomes e 5 (rj) = E^ b g Pb (\y\ < 0.5,po)po/(wR'j 3b ). The charged final multiplicity 
is obtained as 



ay 3 3.6 ay 3 3.6 



1/4 



(10) 



= l^[^RWi{E? b g Pb (\y\<o.5)Y 

= 914 DOl, 
= 1350 DO', 
= 2180 H, 
= 3360 D . 

It is interesting to notice how crucial the lifetime r c of the (Q)GP-phase is for the 
high multiplicities. From Eq. (10) it is observed that due to entropy conservation 

dN ch 



dy 



i ! '* ± C.'C 



since T c is in principle fixed by the theory. Therefore, an increase in the initial parton 
production at r« = 0.1 fm enhances the lower limit of the final multiplicity through 
increasing the lifetime of the plasma. 
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4. CONCLUSIONS 



I have discussed semihard transverse energy and parton production in the mid- 
rapidity unit in central PbPb collisions, and, considered initial conditions for early 
QGP formation at r = 0.1 fm [13]. The main emphasis of the talk has been on the 
consequences of the new HERA parton distributions [15,16] to minijet production at 
the small scales pr ~ 2 GeV. Especially, this lowest order pQCD calculation shows that 
the early QGP is clearly gluon dominated both at RHIC and LHC. For LHC, with the 
new, small- a; enhanced parton distributions it was also shown that rapid thermalization 
of the gluon system may be possible: initially at r = 0.1 fm there are enough gluons 
produced so that e g /n g ~ 2.7T eq as in an ideal gas of massless bosons. However, the 
isotropization of the momenta what was not studied in detail [22]. 

Uncertainties in the pQCD calculation of the initial conditions were pointed out. 
Especially corrections due to shadowing, both in a proton and in a nucleus, were dis- 
cussed, but not yet implemented. For RHIC it seems that the shadowing in a proton may 
be neglected but nuclear shadowing should be taken into account. For LHC, since prob- 
ing smaller x's, one should include all shadowing effects simultaneously [30]. Clearly, 
in this sense, RHIC will be unaffected by the problems and uncertainties in the parton 
distributions at a; < 0.01, whereas LHC will be affected by these complications. How- 
ever, besides problems, the small- ir enhancement causes also larger initial energy and 
particle densities, resulting in longer lifetime for the QGP and larger multiplicities. A 
lower limit for multiplicites in the LHC nuclear collisions was also estimated [13]. When 
the forthcoming colliders are in operation, we will have a truly interesting sequence of 
experiments, AGS,SPS^ RHIC^ LHC to study, and a wide range of multiplicities to 
understand, hopefully in terms of a forming quark-gluon plasma. 
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